Quantum tricriticality and quantum phases in the Lipkin-Meshkov-Glick model 
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We analyse the Lipkin-Meshkov-Glick model near, but not in, the thermodynamic (TV — > oo) limit, 
for arbitrary strengths of the two types of inter-fermion interaction, and at non-zero temperature. 
We use semiclassical methods to obtain the crossover temperature, T*(V,W), between linear and 
nonlinear behavior. It is shown that T* typically vanishes logarithmically as the quantum phase 
transition from the symmetric to the distorted state is approached, except near the quantum tri- 
critical point (V, W) = (0, — e), where it becomes linear. We discuss the relationship of the T* line 
to those observed in conventional condensed matter systems. These calculations are supplemented 
by a direct quantum solution near the line V = 0. This reveals a quantum singlet phase for large 
positive W that is not accessible via semi-classical techniques. 

PACS numbers: 21.60.-n, 21.60.Ev, 71.10.Hf, 71.27. +a 



Introduction. The Lipkin-Meshkov-Glick (LMG) 
model has long been the subject of study in the nu- 
clear physics community. It was introduced in 1959 by 
Fallieros P as a model of monopole oscillations in the 
16 nucleus, and subsequently studied by Volkov [2], be- 
fore springing to prominence with the work of Lipkin, 
Mcshkov, and Glick in 1965 [3H5], who considered it as 
a non-trivial correlated model against which various ap- 
proximation schemes could be tested. 

The Hamiltonian of the LMG model is: 
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where the fcrmionic annihilation operators, er = 

±1 is a spin- like index denoting the nuclear shell in which 
the fermion is, a represents the opposite spin to a, and 
p = 1,2, ...,N P is an auxiliary quantum number dis- 
tinguishing between a large number of degenerate levels 
within each shell. The W interaction represents an ex- 
change of particles between the lower- and higher-energy 
shells, while V represents pair-tunnelling between shells 
(a sort of Josephson term). This Hamiltonian clearly 
conserves the number of particles, i.e. 
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is a good quantum number. It is sometimes stipulated 
that N = N p (the 'half- filled' case), but we shall consider 
all possible values of N. 

The literature on the LMG model is extensive [5] . The 
zero-temperature phase diagram has been obtained in 
the thermodynamic limit [THS] and finite-size corrections 
analysed [10] . Non-zero-temperature properties of the 
model have also been studied [TT], as have properties of 



the zero-temperature entanglement entropy |12) and the 
negativity |13) . Frequently these studies confine them- 
selves to a particular line or region in the two-parameter 
space (V, W) [IH[T5], but some works consider the whole 
plane. It should also be noted that the V = line of the 
model is equivalent to the Dicke model [T|>] , about which 
much is known [I7Hl9| . 

However, there has to our knowledge not yet been a 
study of the finite-iV, non-zero-temperature properties 
of the model over the whole (V, VF)-plane. In this Letter 
we carry out that study, with particular emphasis on the 
nature of the crossover from linear to nonlinear dynam- 
ics that occurs as the quantum phase transitions in the 
(V, VF)-plane are approached at non-zero temperature. 

The large-iVp limit must be taken with care. Since 
the interaction terms in ([I]) cause every level to inter- 
act equally with every other (producing a result ~ Np 
provided that N/N p is finite), we must compensate by 
sending V and W to zero in the following way: 
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with V and W held constant as N p — > oo. Then the 
energy remains proportional to N p , with sub-dominant 
corrections 0(1). However, as we shall see below, these 
sub-dominant pieces are crucial in breaking a degener- 
acy in N p in certain regions of the ground-state phase 
diagram. 

A key observation, made in the original paper [3], is 
that the Hamiltonian may be rewritten in terms of 
the following pseudospin operators: 



'pi) 



J = A. (4) 
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These obey the standard angular momentum commuta- 
tion relations. (Note that we have adopted units in which 
h = 1.) The pseudospin version of the Hamiltonian may 
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easily be shown to be: 



H = e,J z + 



V 



Ji + Ji 



w 

2~K 



(J+J- + J-J+-N), 
(5) 

where N is the particle- number operator defined in ([2]). 
This form of the Hamiltonian is most useful for direct 
quantum treatments; for classical and semi-classical ap- 
proaches it is preferable to rewrite it again using the def- 
initions J± = J x ± iJ y : 
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As well as the particle number, TV, this Hamiltonian also 
clearly conserves the magnitude of the pseudospin: 



Jl + J% + J| = J(J+1), 



(7) 



where J is an integer between and N/2. (We assume 
here and henceforth that N and N p are even.) 

In the remainder of this Letter, we shall summarize 
the existing classical and semiclassical analysis of the 
model, and extract from the latter an explicit form for the 
crossover temperature T*(V, W). Then we shall show, 
by solving the model beyond the semiclassical approxi- 
mation near the V = line, that the ground-state phase 
diagram contains at least one purely quantum phase, not 
captured by the semiclassical analysis. 

Classical phase diagram. Firstly we shall summarize 
the analysis of the Hamiltonian (|6| in the limit J, N — >• oo 
[7j. Note that these limits, while not independent, are 
nonetheless not the same: J —> oo implies N — > oo, but 
not vice versa. In the J — ¥ oo limit, the operators may be 
replaced by classical vectors, and the energy minimized 
straightforwardly. We choose the following parameteri- 
sation: 



•J n 



(1 — cos a) , 



(8) 



(J x ,J y ,J z ) — (J sin 9 cos (f>, J sin 9 sin cj), J cos 9) . (9) 

The maximum possible value of J, denoted J max , is of 
course a function of N: 



J may: 



N 
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< N < N, 



N 



N p < N ^ 2N p . 
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With this parameterisation, the energy as a function of 
a, 9, and 6 is found to be 



hi = (1 — cos a) 
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Minimising this energy, one obtains the ground state 
phase diagram shown in Fig. [l] The three phases shown 




Phase II 



Phase III 



FIG. 1: The ground state phase diagram of the LMG model 
in the classical (N, J — > oo) limit. The bold transition line 
at V — is first-order, while the other two are second-order; 
they meet at a tricritical point when (V, W) = (0, — e). 



in the diagram are characterized by the following behav- 
iors of the pseudospin J: 

Phase I. Full spin, oriented in the negative z direction. 
This corresponds to the parameter values a = it, 9 = tt, 
4> undetermined. This ground state is non-degenerate. 

Phase II. Full spin, canting from the negative z direc- 
tion (near the transition line) to either the negative or 
positive x direction (deep in the phase) . This corresponds 
to the parameter values a = tt, <f> = or tt, and 



W + V 



(12) 



This shows in particular that the transition is second- 
order: there is no jump in the spin's angle of orientation 
as the boundary between phases I and II is crossed. In 
phase II, the ground state is always doubly degenerate. 

Phase III. As phase II, but with <f> — ±ir/2. 

The phase boundary between phases II and III is first- 
order, since while a and 9 are continuous across it, it 
involves a discontinuous jump of the parameter (f>, cor- 
responding to a reorientation of the spin from the x- to 
the y-axis. However, it is a peculiar sort of first-order 
transition, since at the transition all values of <j> become 
degenerate, and hence despite being first-order it does 
have associated soft modes. This emergent U(l) symme- 
try at V — is just the phase of the coherent photon 
field in the superradiant phase of the Dicke model, with 
the quantum tricritical point at (V, W) — (0, —e) corre- 
sponding to the superradiance transition. 

Crossover to non-linear dynamics. We now proceed to 
analyse the finite- N model at non-zero temperature. It 
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is natural to choose a Holstein-Primakoff representation 
[2TT - f23] of the pseudospin, which we define with respect 
to phase I, i.e. with reference to a full spin oriented in 
the negative z-direction: 



J* 
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J_ 
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where the boson operators b and W obey the usual com- 
mutation relations [b,^] = 1, and the commonly used 
linear approximation has been made. Although J is for- 
mally variable, we shall here treat it as fixed at N p /2. 

Substituting this approximation into Q , and applying 
the commutation relations for the 6-operators, we obtain 
that 



H 



eNp 
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V 



+ (e + W)tfb+-({tf) 2 + (b) 2 ). (14) 



It is clear from (14) that the tricritical point (at V = 0) 
corresponds to the point where the boson energy becomes 
negative, signalling an instability which mathematically 
invalidates the linear approximation, and physically cor- 
responds to the superradiance transition. To extend the 
analysis to non-zero V, we must make a Bogolyubov ro- 
tation [24] to eliminate the anomalous terms b 2 and (b>) 2 ; 
this rotation is given by 



b = ,3 cosh 77 + ft sinhry, 
¥ = (3 sinh r] + ft cosh rj, 



(15) 
(16) 



which has been designed to guarantee that [/3,/3^] = 1, 
i.e. that the new operators are also bosons. In terms of 



these, the Hamiltonian ( 14 ) becomes 



H = ((s + W) cosh(2j]) + V smh(2rf))ftp 

+ 1 ((e + W) sinh(2r/) + ^cosh(2?7)) ((3 2 + (ft) 2 ) 
+const. (17) 

The angle 77 that eliminates the anomalous terms is given 
by 



tanh(2r?) = — 



V 



e + W 
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Note that this equation has a solution only in the interval 
|e + W| ^ \V\, i.e. in the area labelled 'phase F in the 
classical analysis above. In the cases in which the Bo- 
goliubov rotation is possible, the resulting Hamiltonian 



H = Eft /3 + const., 
where the boson energy E is given by 
E = ^/(e + W) 2 - V 2 . 



(19) 



(20) 



It is a familiar feature of quantum critical theories 
f3J [55] that the approach to a quantum critical point 



at non-zero temperature is accompanied by a crossover 
from 'renormalized classical' to 'quantum critical' behav- 
ior. A phenomenon of the same sort takes place here: 
the dynamics of the model cross over from being approx- 
imately linear to fully nonlinear as the transition line is 
approached. A simple way to obtain the location of this 
crossover is to ask at what temperature the linear ap- 
proximation (i.e. the condition that (b^b) <C N p ) breaks 
down. The temperature at which this happens is approx- 
imately given by (tfb} T = N p ; inserting the expressions 
for b and b^ in terms of (3 and ft this becomes 
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The value of the thermal average follows directly from 
the Bosc-Einstcin distribution, so that 
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1 



W) 2 - V 2 
which yields the temperature 



= JV„ 
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In this last line we have made several useful definitions: 
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W 
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(24) 

Notice that in the (x, y) co-ordinate system the tricritical 
point is at the origin. 



The first thing to observe about (23) is that it van- 



ishes not at the original phase transition y = x but at 
y = "fx. This is a renormalization of the position of the 
transition line due to quantum fluctuations, similar to 
those discussed in 27J. To examine the behavior of T* 
as this renormalized transition line is approached, we set 
y — jx + i5, with < i5 <C 1. In this limit, provided that 
x ^> 5. we obtain 
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(25) 



as 5 — > + . Hence this second-order transition has an ex- 
tremely narrow quantum critical cone, in contrast with 
the simple power laws typically observed in quantum crit- 
ical theories 25 , 26 . The expression ( 25 ) becomes invalid 
as the tricritical point at (x, y) = (0, 0) is approached. It 
crosses over to a much simpler behavior, which may be 
obtained by setting x — and then taking < y <£L 1: 
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as y — > + . Thus, perhaps surprisingly, power-law be- 
havior is recovered at the tricritical point y — x = 



4 



7" 




FIG. 2: The crossover temperature T* as a function of x and 
y, the rescaled and offset interaction parameters of the LMG 
model. The quantum tricritical point is at the origin; note 
the crossover from logarithmic to linear behavior of T* as 
this point is approached. The line x — corresponds to the 
Dicke model. 



despite being absent along the rest of the critical line 
y = jx. Analysis of the transition from the other phases 
may be carried out in a similar fashion, using the form 
of the Holstein-Primakoff transformation appropriate to 



the phase in question. The full behavior of ( 23 ) is shown 
in Fig. [2] 

Additional quantum phases. To go beyond the semi- 
classical analysis, and thus describe the behavior at tem- 
peratures T > T* , requires in principle a full quantum 
solution of the model. While this is beyond the scope of 
this Letter, it is appropriate to give a brief demonstration 
that the quantum solution of the LMG model contains 
at least one phase that does not appear in the classical 
phase diagram 20 . To expose it, let us slightly rewrite 
the Hamiltonian ([6| with V = as 



H = eJ z 



(27) 



Classically, it is possible to make J 2 — J 2 zero by orienting 
the spin along the positive or negative z-direction; quan- 
tum mechanically, however, this zero is achieved only in 
the singlet state, where J 2 = 0. 

This effect manifests itself at large positive values of 
W/e, where to a first approximation one may neglect 
the e term entirely. The Hamiltonian is then minimized 
by (a) choosing a singlet state for the pseudospin, and 
(b) manufacturing this singlet state from the largest pos- 
sible number of particles; in this case, that number is 
N = 2N P) corresponding to full occupation of all the 
levels in the original LMG model. The energy of this 
singlet state is therefore simply -Egmgict = —W. By com- 
parison, the energy of the phase I state at V = is 
-Efuli-spin — —sNp/2. Hence the transition from phase I 



to the singlet phase occurs when W/e — N p /2. It thus 
recedes to infinite W in the thermodynamic limit, but is 
present at large W for any finite N p . (Further analysis 
shows that there are no intervening phases; the transi- 
tion occurs directly from full- to zero-spin, though via a 
rather interesting quantum critical point.) 

For non-zero V , quantum fluctuations shift the critical 
value of W; to leading non-zero order in perturbation 
theory, the resulting behavior is given by: 



V 2 



The fact that W c increases with increasing V represents 
a pseudo-entropic favouring of the full-spin state, of the 
type recently discussed by Conduit et al. [25] , 

Summary. In this Letter, we have considered the non- 
zero-temperature behavior of the Lipkin-Meshkov-Glick 
model at large but finite N, with specific emphasis on 
the crossover from linear to nonlinear dynamics in the 
vicinity of the model's zero-temperature quantum phase 
transitions. We have obtained the functional form of this 
crossover surface, above which a full quantum treatment 
of the model is necessary. We have also demonstrated 
the importance of a full quantum treatment at large W, 
where a singlet phase emerges for finite values of N p . Fur- 
ther work should be done to explore the properties of the 
region T > T* , particularly near the quantum tricritical 
point. This is especially relevant given the importance 
of tricritical points in the theory of nearly ferromagnetic 
metals [25]. 
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